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Considering the adjacent spin operators SjSj+i and S~ S~ +1 in the 5=1/2 Heisenberg chain, 
we give a determinant representation of their form factors. The dynamical structure factors of 
the respective operators are computed over the whole Brillouin zone in several magnetic fields and 
the resulting signal is analyzed in terms of excitation types. Among other results, we find that 
the SjSj +1 dynamical structure factor carries a large weight of the 4-spinon excitations which are 
distinguishable from the 2-spinon signal because they are located outside the 2-spinon spectrum. 
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I. INTRODUCTION 

One-dimensional (ID) systems are particularly prone 
to quantum criticality as compared to higher-dimensional 
ones. They offer a prolific playground for low- 
temperature phenomenology [TJ [2]. Realizations of ID 
critical quantum liquids are numerous and can take the 
form of e.g. stripes in cuprate high-temperature super- 
conductors, trapped ultracold atomic gases and carbon 
nanotubes. One of the simplest strongly correlated mod- 
els, which remains of great interest due to its underlying 
richness, is the Heisenberg S = 1/2 antiferromagnetic 
chain [3 with Hamiltonian given by 

N 1 - + + - ( z z X \ 

(1) 

where is the number of lattice sites, periodic boundary 
conditions are understood. Several crystals are known 
to realize ID spin chain e.g. KCUF3, SrCuOs, CuPzN 
and have been probed by neutron scattering [4UTU]. For 
a positive coupling constant J > and zero external 
magnetic field ft = the ground state is antiferromag- 
netic without long-range order. Its basic excitations are 
spinons, fractionalized spin excitations that emerge in the 
critical state [11 . The Heisenberg model is integrable 
with an exact solution given by Bethe Ansatz [I2j [13] . 
The properties of its ground state have been extensively 
studied in [14-16 . Moreover, in the case where a spin- 
exchange anisotropy is introduced, one can obtain a com- 
plete phase diagram in terms of the magnetic field and 
the anisotropy parameter [17 . The isotropic Heisenberg 
chain belongs to the universality class of Luttinger liq- 
uids and the asymptotics of its correlation functions can 
be calculated by effective field theory methods. Whereas 
the computation of correlation functions was for a long 
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time a limitation of the Bethe Ansatz method, significant 
progresses have been made in the last decade [HI [19] al- 
lowing the calculation of single spin operator dynamical 
structure factors (DSF) in the 5=1/2 spin chain [20l - 

Eg. 

Going one step further from what has been studied so 
far, we consider the calculation of what we shall call the 
adjacent spin operator DSF 

1 N r°° 
S ab ~ ab (q,uj) = — / dte- iq V- j,> ) +]ut 

i,r=i J -°° 



yields, for eigenstates belonging to the M-subspace: 
\{k} M ) = E (-l) [P1 exp(i Yl v{k Pa ,k Pb )\ 

Pe-KM \ l<a<b<M J 

J2 exp f i^A: Pa j a j |ji,...,j M ) • (3) 

{k} is a set of M quasi- momenta determining the eigen- 
state and the scattering phase (p(k a: kb) is defined by 



2cot^%^ =cot^ 



(S]{t)S b j+1 (t)S^(0)S b jf+1 (0)) . (2) of M localized down-flipped spins as 



9 cot ^ . We introduced the state 



where a, b = z,—,+. Resuming the study of the spin- 
exchange DSF in [29], our interest lies in the computation 
of specific components of the latter correlation function. 
In addition to the resonant inelastic x-ray scattering re- 
sponse function studied in this previous paper, the adja- 
cent spin operator form factors (FF) are related to the 
study of transport and spin currents along the spin chain 
[3QH34] and to the fidelity in the Heisenberg spin chain 
[23 [36]. 

This paper presents a non-perturbative calculation of 
the adjacent spin operator DSF for (a, b) = (z, z), ( — , — ) 
(in (J2|) in the antiferromagnetic S = 1/2 Heisenberg 
chain. 

Low-energy effective field theories have already been 
used to predict the asymptotics of the adjacent spin op- 
erator correlation function in a few cases [33 [38] . 

The contents of this publication are the following. In 
section [II] we present the exact eigenstates of the system 
provided by the Bethe Ansatz. Using a determinant rep- 
resentation, the FFs of the operators SjSj +1 and S~ S~ +1 
are formulated in terms of ei genstates rapidities in section 
III We show in section IV and |v| how the S h+ (g,co>) 
and S Az (q,uj) DSF are constructed using the FF expres- 
sions. In section \VJ\ a catalog of the excitations is given 
by describing the quasiparticles which occur in presence 
or absence of magnetic field. In section |VII[ we evalu- 
ate numerically the expression for the two DSFs over the 
entire Brillouin zone. The contribution of each family of 
excitations is shortly discussed. The results are summa- 
rized in the conclusion in section lYUH 



II. SETUP 



(4) 



It is convenient to express the quasi-momenta in terms 
of rapidities A, 



exp(ifci) 



Xi + i/2 
A< - i/2 ' 



(5) 



Imposing periodic boundary conditions leads to the 
Bethe equations for the rapidities 

7T 1 M 

arctan (2A$) = —Ii + — ^ arctan(A^ — A&), Mi . (6) 



N N 

k=l 

Taking N to be even, the quantum numbers 7i, . . . , Im 
are a set of distinct integers for odd M and half-integers 
for even M. Each of these sets of quantum numbers 
specifies a set of rapidities and vice versa. It can be 
shown that the set of quantum numbers for the ground 
state is {I® = i— M 2 ^ 1 }, i = 1, . . . , M. Once the rapidities 
of an eigenstate are obtained, the energy and momentum 
of the state follow straightforwardly 



M 



E 



1/2 



N 

2 K--M) 



p = ttM - ^ ^2 Ii (mod 2tt). 



M 



N 



(7) 



The Bethe equations represented in (|6| provide a sys- 
tematic way of constructing eigenstates and the solutions 
allow one to obtain dynamical correlation functions and 
thermodynamic quantities with high precision. 



A. Bethe Ansatz 



B. String solutions 



Thanks to the integrability of the Heisenberg spin 
chain, the exact eigenfunctions of (TJ can be obtained 
via the Bethe Ansatz [12] . Because the Hamiltonian con- 
serves the total magnetization of the chain, the Hilbert 
space separates into subspaces of fixed magnetization 
M = N/2 — (^2^ =1 Sj). Defining a reference state with all 
spins pointing upwards |0) = ®^Lilt)> ^ ne Bethe Ansatz 



The solutions of §6y are not restricted to real rapidi- 
ties only. As Bethe mentioned in his seminal work [12], 
self-conjugate complex solutions, which can be inter- 
preted as bound states of down spins, also satisfy the 
Bethe equations. Bethe [12] made the conjecture, com- 
pleted later by Takahashi et al. [39 , that such com- 
plex rapidities form string structures which are A^' a = 
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X J a + |(nj + 1 — 2a) + i(%' a with a = 1, . . . , rij, rij being 
the string length and the string deviation 5a J,a . Here a 
is an index that runs from 1 to Mj where Mj is the num- 



ber of rij strings (strings of length rij) and j 



where N s is the total number of different possible lengths. 
The string hypothesis is unfortunately not always cor- 
rect, although it can be shown that generally the string 
deviations vanish exponentially with system size. The 
correct number of complex solutions is moreover smaller 
than the one constructed from the string hypothesis and 
other kinds of complex solutions appear. The propor- 
tion is of order 1/y/N [40]. More extended discussions 
about string structures can be found in [41]. The devi- 
ated string solutions turn out to be unimportant for the 
correlation functions presented in this paper. 

In the presence of string solutions with vanishing de- 
viations the Bethe equations ([6| become undetermined. 
The remedy is to rewrite the Bethe equations only in 
terms of the real centers A« J of the strings to obtain the 
Bethe-Takahashi equations [39] : 

N s M k = l 
k=l (5 



The norm of a state which contains one or two infinite 
rapidities can easily be calculated by commutation of spin 
operators and are 



7V({A,oc} M ) 
N({\, oo,oo} M ) 



N - 2M + 2 
N 



7V({A} M -i) , (10) 



2 (TV — 2M + 4) (TV — 2M + 3) 



N 2 



• N({X} M -2) 



(11) 



III. FORM FACTOR DETERMINANT 
REPRESENTATIONS 

The natural language for expressing normalized eigen- 
states and form factors (FF) is the algebraic Bethe 
Ansatz [I9j [2TJ [22] . In the algebraic Bethe Ansatz for- 
malism, one introduces the monodromy matrix T(A) that 
acts on the space C 2 <8>Hn where Hat is the Hilbert space 
of the TV-sites chain. The operator is represented in the 
auxiliary space C 2 as 



_ (A{\) B{\)\ 
T ^ X) ~ \C{\) D(X) 



(12) 



where 

nj (X) = 2atan(2A/n i ) 

B Uj , nfc (A) = (1 - Sj,k)0\ nj - nk \ + 2(X) + 20 nj - nk - 2 (X) + 

+ ... + 20 nj + nk - 2 (\) + nj+nh (\) (9) 

and I£ are the set of quantum numbers corresponding to 
strings of length rij . 



where A is called the spectral parameter. The 
monodromy matrix is constructed as a product of 
L— operators 



T(A)=L iY (A;£ iV )...MA;£i) 



(13) 



with Lj(\; £j) = Roj(X — ) ' wnere refers to the auxil- 
iary space and j to the j— th site of the spin chain. The 
R— matrix R^ must be a solution of the Yang-Baxter 
equation. The £j are inhomogeneity parameters and the 
homogeneous spin chain corresponds to setting all 
C. Spin of the eigenstates and infinite rapidity ~> V 2 - However, it is more convenient to take this 

limit at the end of the computation, as we will do in 
this work. One can now define the transfer matrix as 
T(X) = A(X) + D(X), from which the Hamiltonian can 
be recovered 



In the isotropic spin chain the total spin operator 
(SjSj) commutes with the Hamiltonian and in con- 
sequence its eigenvalues S tot are conserved. More- 
over, the eigenstates created with the Bethe Ansatz are 
highest-weight with respect to the global su(2) algebra 
i.e. the states constructed with M rapidities have the 
spin eigenvalues = S tot = f - M g2]. In or- 

der to access eigenstates with S l ° l < S tot , one must 
act on a Bethe state with the global spin lowering op- 
erator S~ =0 — -^^^2^=iS~. These descendant states 
are actually also Bethe solutions but they include in- 
finite rapidities. Indeed the Bethe equations allow ra- 
pidities to go to infinity and from eq. [3] one notices 
that a state with M rapidities of which M' tend to 
infinity can be written | { Ai , . . . , Xm-M' , oo, . . . , oo}) = 



H 



2 dX 



lnT(A) 



(14) 



X=i/2 



Hence the transfer matrix and the Hamiltonian have com- 
mon eigenfunctions and can be constructed as 



B(Ai)...B(A M )|0) 



and the dual state as 



<0|C(A 1 )...C(A M ) 



(15) 



(16) 



[S q=0 ) lvl |{Ai, . . . , Am-M'})- This eigenstate is no 
longer highest-weight but has S tot 
St ot = % - M. 



f - M + M' and 



where the rapidities Xj have to satisfy the Bethe equa- 
tions ©. The norm of eigenstates can now be expressed 
as a determinant [42-44]. Another very import result is 
a determinant expression for the scalar product [44] 

(0|C(/i!) . . . C(/z M )fl(Ai) • • • B(X M )\0) (17) 
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where the fij can be arbitrary and the Xj have to satisfy 
the Bethe equations or vice versa. A last essential ingre- 
dient is the inverse mapping that expresses the local spin 
operators aj,<7^,<7~ in terms of the non-local operators 
A(X),B(X),C(X),D(X) EI] 

j-l N 

<r! = -2l[[A + D]{t i )D(t j ) II [A + D](£ k ) + 1, 

i=l k=j+l 
j-l N 

^ =n^ +D ] n i a + d ](^)i 

i=l k=j+l 
J-l N 

<=n[ A+D H6)c(ei) n [A+Duth). as) 

k=j+i 



i=i 



j 



Using these expressions and the scalar product (17), we 



expressed in terms of determinants the adjacent spin op- 
erator form factors which are defined as 



(GS\S«S b j+1 \a 



(19) 



with \GS), \a) respectively the groundstate and an eigen- 
state of the system. The form factors involving a single 
spin operator were first derived in [21], the general case 
was treated in [22]. These expressions are not always 
suitable for a numerical evaluation, therefore we re-derive 
below the form factors of interest for this paper. 



A. SjSj +1 



With ({A} | and \{fi}) two eigenstates of the system satisfying the Bethe equations, the form factor of the operator 



qz qz 
3 ^'+1 



4^ + 1 



reads |({A}|cr^<j^ +1 |{/i}). Using (18) and the properties of the transfer matrix, the product of 



two adjacent o z , operators is 



3-1 N 

i=l i=j+2 
+<7?+<T? +1 -l. (20) 

The form factor for a single cr| operator has been considered in [ 2TJ [25] and the remaining non-trivial part of cr|<j| +1 
form factor is the term: Yli=i i A + D] (£i)D(£j)D(£j+i) Ylk=j+2 i A + D\ Using the expression for the action of 
the D(X) operator on a general state [19] twice, we can write the form factor as a double summation of Slavnov 
determinants 



with 



n=l i^n ^i) m=l,m^n i^m + 1 ^i) 



A2 



(v) 



4>{^m -fJ-n + V) HfJ-m, Ul>k H X k ~ A;) Ul<k HVk ~ Mi) 

det H({Xj}, {/J,ijt m ,„, £j , gj+i }) 



i+U 



(21) 



H ab ({X}, {/x}) = 0(A ^ b) I n <P(X k -to + v)- d(ji b ) J] ^(Afe -Hb-V) 

d(X) = f\b(\,£i), b(fx, A) = ^'x + r,) ' ^ 

The value of r] and the definition of 0(A) depend on the anisotropy parameter of the spin chain [2Tj . We keep their 
general expressions for the rest of this section in order to cover any anisotropy. However, in the isotropic case they 
take the values 



0(A) = A, r] = i 



(23) 



5 



1. Homogeneous limit 

The homogeneous spin chain corresponds to the case where & — »• 77/2 Mi. This limit should be taken with care in 
order to obtain a finite expression for the ratio deti^({A}, {/i^ m , n , £7, £j+i})/<K£j+i — We use l'Hopital's rule to 
have a well-defined homogeneous limit and the matrix elements become 



{{X}\D( V /2)D( V /2)\{^) = 



M M 



J2 An Yl detpr ' 



rn=l 



(24) 



Here we defined 



M 



A n 

nn 
B ab 

c a 



rrinm 
*ab 



G 



ab 



d(/i n )0(/i n - 77/2) Y[ <t>(l*>i - Vn-rj) 

i=l 

(1 - (5 6>n )d(/i 6 )0(/i6 + 77/2) Yl <t>(l*>i ~ Vb-v) 

2 (A a -r7/2)0 2 (A a + 77/2) 
b = m 



0(77) 



0(A a - /7/2)0(A a + 77/2) 



4>{\ a —ii h ) 



(rifc#a ^ - + ^7) - d(/i 6 ) ri/c/a 0(Afe ~ »b ~ Tj)) , h + 



n 

b = n 



(25) 



Using Laplace's determinant formula (see chapter 9 appendix in [19 j) and the fact that B n is a rank one matrix, we 
write the summation over determinants as a single determinant: 



A I 



^detF nm = det(G n + B n )-detG n 



(26) 



and the complete form factor is then 

<{A}k>| +1 |M) = <{A}|cr||M) + <{A}|a| +1 |{ M })-({A}|M) 
, A <Pj-i(W) n^ 2 (A i + r?/2)r 2 (Mi + ^/2) 



ft I 



vj-i(M) n i>fc 0(A fe - \j) n i<fc - Mi) ^ 



^A n (det(G" + B n )-detG n ) . (27) 



Here we introduced 



and if {A} is an eigenstate, then we can write — e _ljP ^ A > with P{\] the total momentum of the state. 



(28) 



2. Fourier transform 



For the further computation of the dynamical structure factor, the Fourier transform of the form factor (27) is 
necessary. With (1/y/N) Ef=i ({A}|5J5J +1 |{/x}) = (1/y/N) E P e ip ({A}|S|_ p S*|{/i}), we give here an explicit 
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expression for the norm squared of theses matrix elements between normalized states 



2tt 



E p e ip ({A}|^_ p ^|{ M }) 



N 



N({\})N({»}) 



Fi 

F 2 



H ah ({\},{»}) 
P a b({A},{ M }) 



N 
16 



n 



<t>{\ + J7/2) 



<t><jH + 7?/2) 



n a >^(^(Aa-A /3 ) + ^( r? ))|$( {A })| 

|( e i(9A-9 M ) + 1) i? x +e i(9A-9 fl ) j P 2 _ ({A}|{^})F 3 | 2 

det(F({A},M)-2P({A},{ M })) 

- 4 n£^|E^(^ + ^)- d et G ») 



j>k lLiY 



<t>iji) 



0(A a - Mb) 



(Ai + 77/2)0-2(^ + 77/2) 



cj ) (X a -7 1 /2)cj ) (X a + V /2) 



ky^a 



(29) 



with iV({A}) the norm of the state [42, 43, 45 including the Gaudin determinant |$({A})|. Recall that 0(A) = A and 
77 = i for the isotropic spin chain. 

In the case that the state ({A}| contains string configurations, columns of the matrix become equal to leading 
order in the string deviations 5 and one should use the reduced Gaudin matrix |\I/( r )({A})| instead (see j46], [25]). 
The other factors in (29) are well-defined if only the state ({A}| contains strings. However, when the state 
contains string structures, the entries of the matrices G n , B n , H, P yield undetermined factors. In case of a single spin 
form factor, the matrices H and P can be regularized (see [25] and the discussion below ([33]) ) . Unfortunately these 
methods do not apply for the form factor under consideration here. But since the operator SjSj +1 is Hermit ian and 
the ground state is made of real rapidities, one can always choose ({A}| to be the state with strings, and obtain well 
defined results for the form factor. 



B. Sj S j+1 



The second FF of interest concerns the operator Sj = Oj cfj +1 . Again, using the result ( [18] ) for the local a- 



operator and simplifying the product of transfer matrices, the form factor reads 

i-i at 
{{X}uWjaJ +1 |{/i} M - 2 ) = {{\} M \X[(A + D)^)B(i 3 )B(i 3+1 ) \{ (A + D)(£k) \{^}m-2) (30) 

2=1 k=j+2 



with the use of the determinant representation and taking the homogeneous limit similarly to (25), the form factor 
becomes 



^-i({A}) Elf <KA, + 77/2) 2 detH— ({A} M , {mJm- 2 ) 



({A}m|<^ (7j + 1 |{/i}M-2) 



with the H matrix 



<p(X a -fi b ) 



vj-im) nf _2 + v/2) 2 n£* ^ - \j) Y\f- k 2 - h) 

(a(»b) Uk^a 0(A* - lib + 77) - d(n b ) Uk^a ^(Afe - »b - 77)) , b^M-l,M 



(31) 



0(A a -r 7 /2)(/>(A a +r ? /2)' 
0(r7)0(2A a ) 
I 2 (A a -77/2)02 (Aa+r7/2 ) 5 



b = M- 1 

6 = M 



(32) 
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1. Reduction of determinant for string solutions 



If one considers states which include string structures, following the same procedure as [25 j, the determinant 

representation ( [32] ) has to be reduced in order to have well-defined matrix elements. Hereafter, we use two different 
ways of indexing the rapidities. The notation with one index /i a (a = 1, . . . , M) allows one to identify a row or a column 
of the matrix under consideration. On the other hand, the notation including three parameters \j^ a (k — 1, . . . ,N 3 , 
a — 1, . . . , Mk and a = 1, . . . , rik) provides a clear identification of rapidities which belong to the same string. The 
reduced H ( r ) matrix reads then 



H 



ESS: [d-4 



Ji tt> FiFi+i 



4>(ji)K a ,i-iJi-i, 



Si,o)L a ,i - (8i, n + S ifi ) K aA ] , b:(k,a,l),b^M-l,M 
b : (k, a, i), i = 2, . . . , n k , bj^M — \,M 



M 



K„A = 



02(A o -r ? /2)</»2(A o + r ? /2) 
1 



k,i+l^ 



L a ,i — 



d\ a 



(33) 



2. Fourier transform 



In order to express the dynamical structure factor, one needs the Fourier transform of the form factor (32) 
(l/VW)Ef=ie-^({A}|5757 +1 |{/i}) = (l/VN)Z p e ip ({\}\S-_ p S-\{^}). We express here the norm squared of 
these matrix elements between normalized states: 



2tt 



E p e ip ({\}\S-_ p S-\{fi}) 



N 



iV({A})iV(M) 



<t>W2 + Mi) 2 



5{q + q» -q\) 



Ua>p^HK-^) + ^mM{M)\ 

,,2 



F 



(34) 



with 



detH— ({A},{ M }) 



A i)n^ fe 2 0(Mfe "Mi) 



(35) 



and with iV({A}) the norm of the state [42, 43, 45 including the Gaudin matrix determinant |3>({A})|. We recall that 
(j)(X) = A and r] = i for the isotropic spin chain. If the state {/i} contains strings, this results has to be combined with 
the reduced determinant explicitly given in Q and the reduced Gaudin matrix $ (r) (M) (06], [25]). 



IV. S Z jS z 3+1 DYNAMICAL STRUCTURE FACTOR IN THE HEISENBERG SPIN CHAIN 



In order to have the response of the system to an excitation of momentum q and energy cj, we express the dynamical 
structure factor ([2| in Lehmann representation with the Fourier transform form factor: 



N 



(36) 



where (GS\Sq_ p Sp\a) is the form factor between the normalized ground state and the whole set of normalized 
eigenstates and u a = E a — Eqs is the relative energy. We will explain hereafter that in order to cover all the 
contributing eigenstates |a), we express the results in terms of the SjSj +1 and S^Sj~ +1 form factors. 
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FIG. 1: Sr(q,u) DSF for a N = 400 XXX spin chain and (a) M = 50, (b) 100, (c) 150 and (d) 200 



A. Spin sectors 



As the Sg pSp conserves the spin S z and if we consider the form factor with a ground state with M reversed spins, 

the matrix elements are only non-zero for eigenstate with the same magnetization. To access these states in the different 
sectors but with S z ot = N/2 — M, we include the eigenstates that contain one and two infinite rapidities, and with 



the contributions of all sectors, the DSF reads: S 4z (q, u) = Em'=o E a , E P ^ v (GS M \ S z q _ p S z p (S^) 



\(M-M') 



OLM') 



with (GSm\ the normalized ground state of magnetization S° = y — M and \cxm') a normalized eigenstate including 
M' rapidities. As the ground state is also highest-weight, Ei S^~\GSm) — an d by commutation we identify the only 
three contributing sectors 

E|E eip < gs m\ sz a-v s p i«) 2 = E I E eip ( gs m\ s I-p s p m 



- jv_ 2M+2 E Ie^^mIv^+^-a-i^m-!) 

OCM-l P 

2 £ |E e "'< GS «iw;i a «-3) 



(N - 2M + 3)(N - 2M + 4) 



(37) 



OCM-1 P 



where the first and last Fourier transform form factors given in (29) and (34). In the case of zero magnetic field, 



(M = N/2), the contribution of the sector M — 1 is identically zero. Indeed, following the Wigner-Eckart theorem [47] . 
the matrix elements can be expressed proportionally to a Clebsch-Gordan coefficient. We decompose the indexing 
of the Bethe states \a) into a state parameter, total spin and spin-z number: S tot , Sl ot ) and we rewrite the form 
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factor as: 



(GS\ S*S* +1 \a) cx J2 ( G5 .°.°l S j |/M'°',4 ot > (/?,5 to *,4°*|^ +1 |a,S tot ,0> 

^ gtot gtot 



oc 



^ (0, 1; 0, 0\s to \ sl ot ) (s to \ 1; 4 ot , 0|S tot , 0) 



cx (0,l;0,0|l,0><l,l;0,0|5* ot ,0> (38) 
with (ji, j'2; mi, J, M) the Clebsch-Gordan coefficient. Then the only non-zero contribution sectors are S tot = 



N 



2 - M, f - M + 2. 



One notices from ( |37| ) that for large spin chain, in a finite magnetic field, the contributions of the sectors M — 1 
and M — 2 are in order of and , so negligible. In summary, to evaluate the dynamical structure factor S 4z (q,uj) 
one needs the S*_ p S* form factor in sector 5' tot = y — M and for a magnetization close to zero (TV — 2M ~ 1), one 



must add the contributions of S q _ p S p form factor in sector 5*°* = 2. 



B. Sum rules 

The sum rules are of particular importance for the calculation since they provide checks independent of the algebraic 
Bethe Ansatz. The integrated intensities allow one to control the efficiency of the algorithm by summing all the 
computed form factors. 



1. SjSj+i integrated intensity 



N ^ 2ir J ' 



duS iz { q ,u) = ^ EE e " i9(i_/) < GS I s i s Ui s r s r+i l GS > 



This result is valid for every value of the magnetic field and for any anisotropy. 
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(39) 



2. SjSj+i first frequency moment 



J II <i <if 



M 3,3 



N2 



^ e -iqV-?) (GS| 



EQX QX j_ qU QU qz qz 



QZ QZ 



|GS) 



J 

iV 



£ (GS| J (S;Sj +1 + SjSjVi) " cos(2 9 )5;_ 1 (S;SJ +1 + SjSjVi) S; +2 |GS) . (40) 



For h = and in the isotropic spin chain, we use the expectation values calculated in [48] to evaluate the first frequency 
moment that reads then 



J 



(41) 



V. 5757 +1 DYNAMICAL STRUCTURE FACTOR IN THE HEISENBERG SPIN CHAIN 

For two adjacent local spin raising operators, using the result ( |34[ ), the dynamical structure factor reads 

S- ++ (q,cu) = ^y\y^(GS\S-S-\a) 2 5{uj-uj a ) (42) 
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FIG. 2: S— ++ (q,u) DSF inaiV = 400 XXX spin chain and (a) M = 50, (b) 100, (c) 150 and (d) 200 



where (GS\S~_ p S~ \a) is the form factor between the normalized ground state and the whole set of normalized eigen- 
states and uj a = E a — Eqs is the relative energy. The operator S^L p S™ flipping down two spins and by conservation 
of the z-spin value, no sector of M' > M — 2 is contributing to the dynamical structure factor. Moreover following the 
same procedure as (37) by acting with an S~ =0 operator, it is straightforward to show that the contributions of higher 



spin sector for M' < M — 2 are all null and that if the ground state belongs to the Mth sector, the only non-zero 
form factor are with eigenstate of sector M — 2, |qm-2)- 



A. Sum rules 



In this subsection, the integrated intensity and the first frequency moment of the Sj operator are calculated. 



1- Sj Sj +1 integrated intensity 
Summing all the dynamical structure factor over all frequency and momenta gives the intergrated intensity: 



N ■ 

3 

4 ~ 1 2 ~ ~N ) + \N~dA + 4 



(43) 



where we used the Hellmann-Feynman theorem for the evaluation of ^2- (GS| SjSj +1 |GS). The result is valid for 
h > and for the isotropic as well as for the anisotropic spin chain. 
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2. S j S j+1 first frequency moment at zero magnetic field 



The results hereafter are restricted to the case of zero magnetic field since we use several rotational symmetries in 
order to keep the result concise. 



2tt/ 



27V 



^e-M-j') (G S| \[H,S-Sr +1 ] ,S+S+ +1 ] |GS) 



2N 



£-4co S (2<z) <5|5| +1 5| +2 5| +3 ) + 2co S (2q) (S+ S~ +1 S* +2 S* +3 ) 



- (6cos(2 9 ) +4 008(9)) (SfS* +1 S* +2 Sr +3 ) + (4co S (q) - 2cos(2g)) (S+ S* +1 Sr +2 S* +3 ) 

+ (2-2cos( g ))(5^ +2 ) + (4 + 2cos( g ))(S?^ +1 ) • (44) 

We denote here (O) the expectation value of the operator O at zero temperature. This result is only valid for isotropic 
spin chain. 

If we use the results of [48] to evaluate the correlators, the first frequency moment becomes 



J dujujS—(q,uj) = 21n(2)-^C(3)-^+cos(c7) 



C(3)--ln(2) 



cos(2g) 



ln(2) flOC(5) - - - 5C(3) 



?«»> + «.>(H«*> 



1 

20 



(45) 



VI. IDENTIFICATION OF THE EXCITATIONS IN THE ISOTROPIC SPIN CHAIN 
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FIG. 3: Combinations of quantum numbers corresponding to the ground state and several excited states (see text for description) 
in the N = 16 spin chain. 

We give in this section a description of the excitations in terms of quasiparticles that occur in the Heisenberg spin 
chain. This allows one to identify the contributions to each correlations. We divide the description of the excited 
states into two parts: zero and non-zero magnetic field, since the polarization modifies drastically the combinatorics 
of the states above the ground state. The identification of the states is made with the use of the set {I 3 a } (see ([6|,([8])) 
and the excitations of the spin chain are constructed from all the possible combinations of quantum numbers within 
the boundaries 39 



\Pa\<Pcc 



1 



N s 



N + 1 - M k ( 2 min(n j , n fc ) - S jik ) 



k=l 



(46) 



Excitations at h — 



In the case of zero polarization, the ground state quantum numbers occupy all the possible vacancies for finite real 
rapidities, i.e. {/* = a — M 2 ^ 1 }, a = 1, . . . , M. The creation of excitations requires then to either take a rapidity to 
infinity, create a string or to remove rapidities. 
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1. String excitations 



The creation of an n string configuration allows vacancies in the quantum numbers for real rapidities. The string 
state so created is non-dispersive since there is only one accessible quantum number, but the holes created in the set 
denoted "spinon" quasiparticles [IT], have a continuum spectrum of excitations. The presence of the complex 
rapidities imposes via the Bethe equations the number of spinons created and for each n string, there are 2n — 2 
spinons in the excited state. As notation, we use, e.g. 2spls 2 for a state with 2 spinons and one 2 string. 



2. Spin raising and infinite rapidities 



There are two other processes which create spinons by removing quantum numbers from the ground state. 

First, if one acts with a local spin raising operator, i.e.: Sj~, the spin chain will change sector from M = N/2 to 
M' = N/2 — 1 in an excited state with S tot = Sl ot = 1. This is equivalent to removing one number from the set {/^}. 
The quantum numbers go then from integers to half integers or vice versa, and consequently each spin raising creates 
two spinons. We denote this excitation with 2sp, 4sp, ... for two, four and higher number of spinons. 

Second, the Bethe equations allow one to send rapidities to infinity and therefore to remove the corresponding 
quantum numbers. Similarly to the spin raising operation, this creates a pair of spinons for each rapidity sent to 
infinity but on the contrary of the spin lowering operation, the excited has 5^ ot = but S tot = 1. The infinite rapidity 
is actually equivalent to a global spin raising operation: | {A, oo}m) = Y^j $J \{^}m-i) where |{A,oo}m) belongs 
then to a higher spin sector but with the s ame ff * ot . These states are described in detail in section II C and they 
contribute to other spin sectors (see section IV A). The notation we use is 2spoo, 4sp2oo, ... for excitations of one, 
two, etc. number of infinite rapidities. 

For the three different kinds of excitations described above, the quasiparticles created, spinons, have always the 
same dispersion relation and therefore the spectra are identical even though the states are actually different. 



B. Excitations at h > 



In presence of a magnetic field, the ground state of the spin chain is polarized, 5^ ot > 0, and therefore there are 
N — 2M vacancies in the quantum numbers. Similarly to the h = case, the lowest energy state forms a Fermi like sea 
in the quantum numbers: {I* = a — ^^}, ol = 1, . . . , M, although there are vacancies between the highest /lowest 
quantum numbers and the boundaries: M ~ 1 +1 < 1^. We discuss here three excitations that can occur above the 
ground state. 



1. Particle-hole 



Due to the presence of vacancies, a quantum numbers from the ground state can be taken to higher value. One 
creates this way a pair of quasiparticles: a hole left in the Fermi sea and a particle evolving above the Fermi surface. 
We write in this article these excited states by ipih for i particle- hole pairs. 



2. Spin raising 



Similarly to the h = case, acting with a spin raising operator Sj~ on the spin chain, removes one rapidity and 
creates a hole in the Fermi sea. The change in qantum numbers parity change also the number of vacancies above 
the Fermi sea but we only count as quasiparticle the holes below the Fermi surface. Therefore, applying i spin raising 
operators on the ground state creates an excited state with i holes belonging to the sector M' = M — i and which 
we denote ih. In finite magnetic field, the creation of infinite rapidities is still admissible by the Bethe equations and 
provides similar hole quasiparticles but as discussed in IV A the form factors of such states decay like the inverse of 
the spin chain length. 
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3. String excitations 

Excited states can be created by forming a string configuration in the rapidities. In contrast to h = 0, a n string 
leaves n holes in Fermi sea of rapidities. Moreover, the boundary is of O(N) and then much bigger than at h = 0. 
Therefore the string states have dynamics that contribute to the spectrum. The notation we use for a i string state 
is ihls\ 



VII. DSF EVALUATIONS 

In this section we discuss the evaluation of the dynamical structure factors (DSF) S 4z (q,u) and S h+ (g,cj) 
computed following the schemes IV and [V] and the contributions of the different kinds of excitations are analyzed in 



the principal cases. 

Considering a spin chain of length N = 400 and with fillings of M = 200, 150, 100, 50, the value of the DSF is 
computed by summing over the eigenstate contributions. The summation over the intermediate states is performed to 
obtain quantitative results over the Brillouin zone and over an energy range that covers all the significant weight. This 
is done using the ABACUS algorithm [49] which sums intermediate state contributions in a close to optimal order. 
From a given set of quantum numbers, we solve the Bethe equations a nd calculate the energy and momentum 
([7]). The value of the form factor is then computed with the formula ( 29|34 ) . 



The saturation of the integrated intensity is checked along with the first frequency moment in order to control the 
completion of the computation. From a general point of view, the saturation could be improved by either lengthening 
the time o f computation or improving the scan efficiency. To obtain smooth curves in frequency cj, the delta function 
in ( 36|42 ) is broadened in energy to a scale commensurate with the level spacing. 



■ ++ (q,oj) 



(%) 


M = 50 


M 100 


M 150 


M 200 


2h 


98.1 


88.3 


64.3 


4sp 90.8 


lp3h 


1.9 


11.2 


29.8 


6spls 2 1.4 


2p4h 


0.01 


0.4 


2.4 






" ++ 99.9 


99.8 


96.6 


92.2 



TABLE I: Contributions of excitation types to the the S — ++ (<?,cj) DSF for N = 400 and M = 50, 100, 150 and 200. The last 
row gives the saturation of the integrated intensity . 

The S h+ (g,cj) evaluation is shown in figure [2] for the four different magnetizations. One observes that the shape 
is similar to the single spin raising S ^~(q,uj) DSF~|50 but with a shift of tt in momentum because the two spin raising 
operator removes two quantum numbers to the ground state instead of one. There is also a significant difference in the 
distribution of the signal. For the two spin raising DSF, the weight is more homogeneously spread from low to high 
energy although for a single spin raising, it is mainly located at the lower boundary of the spectrum. S h+ (g,cj) 
is the propagator of a pair of neighboring spin up and as a consequence, it vanishes when the polarization of the 
spin chain tends to saturation (M = 0). In table [ij one can observe also that as the magnetic field increases, the 
number of excitations with particle- hole pairs created by the double spin raising decreases. In figure |4j we split the 
total signal of S ^ + (q,uj) into the three main contributing types of excitations. The signal becomes weaker as the 
number of particles in the excited state increases and it is also noticeable that if the 2h type carries weight at low 
energy, the signal of the lp3h and 2p4h vanishes at zero energy and is mainly significant at higher energy. Besides, 
the DSF at h = (M = 200) shows a signal which is mainly composed of 4-spinon excitations (table [i]). We divide 
the S h+ (g,cj) DSF signal into the 4-spinon and 6-spinon excitations in figure [ij In addition to the fact that the 
6spls 2 contribution is globally two orders of magnitude lower than for the 4sp, we notice that the 6spls 2 weight is 
essentially inexistent between q = ir/2 and q = 37r/2. 
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(a) 2h 



(b) lp2h 



(c) 2p2h 



0.018 

0.016 

0.014 

0.012 

0.01 

0.008 

0.006 

0.004 

0.002 







(d) 4sp 



(e) 6spls 2 



FIG. 4: Contribut ions by type of excitation to the S DSF in a N = 400 spin chain for M = 150 in (a) (b) and (c] 

or for M = 200 in (d) and (e) The graphs are labeled by the type of excitation that are represented. 



B. S 4z (q,uj) 

The S Az (q, to) DSF is here numerically evaluated and the signal weight is analyzed and attributed to the different 
kinds of excitations. In order to emphasize the dynamic properties, we consider only the connected DSF defined as 
S* z (q,uj) = S 4z (q,uj), V(g,w) ^ (0,0) and Sf{q = 0,w = 0) = 0. 

S 4 c z (q,to) 



(%) 


M = 50 


M 100 


M 150 


M 200 


lplh 


93.2 


86.4 


64.8 


2spls 2 


57.4 


2p2h 


5.3 


7.0 


14.4 


4sp2oo 


38.0 


3p3h 


0.1 


0.1 


0.2 


6sp2ools 2 


0.6 


2hls 2 


0.1 


1.2 


9.1 


4sp2s 2 


0.003 




98.68 


94.69 


89.01 


96.0 



TABLE II: Percentages of contributions of excitation types to the the Si z (q,u) DSF for N 
200. The last row gives the saturation of the integrated intensity. 



400 and M = 50, 100, 150 and 



As pictured in figure [I] in the presence of a magnetic field the signal is very similar to the single spin S zz {q,uo) 
DSF [25j|50] and as shown in table [n| the contribution of the 2 string states are of the same order. From the table Il| 
one notices also that the percentage of the 1 and 2 particle-hole excitations to the sum rule is approximately stable 
as the magnetic field decreases. From the analysis of the compu tation resu lts, we ca n identify where the different 
kinds of excited eigenstates supply weight and we show in figures 5(a)| |5(b)| and 5(c) the contributions of each type 
of excitation. For the case M = 150, we split the S^ z (q 1 uj) DSF into the three main weight carrying types. The lplh 
states contribute mainly at low energy whereas the 2p2h states carry their highest signal in the upper part of the 
DSF map. In addition, the weight of eigenstates with a single 2 string structure is located around q = tt and uj = 2 
and this category of excitation is gapped in energy. Although the string is accompanied by 2 holes, the signal clearly 
shows the presence of this bound state which is difficult to observe [5T] . 

The signal of the S^ z (q^uj) DSF at zero magnetic field takes the original form shown in figure [l] As shown in 
table [nj most of the weight carried by the spin singlet eigenstates (Stot = 0) corresponds to 2spls 2 excitations. By 
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definition, these excited states are located inside the 2-spinon spectrum, i.e.: ||sin(g)| < cj2-spinon < 7rsin(g/2) [TT] . 
Almost all the rest of the weight is carried by 4-spinons states (6 and higher spinons states are negligible) and falls 
within 4-spinon spectrum: ^|sin(g)| < co>4_ S pinon < max{27r sin(g/4), 2tt sin(g/4 + 7r/2)}. We see on table[ll]that among 
the 4-spinon states, the spin quintuplet eigenstates (with two infinite rapidities) contribute the most to the signal 
and si nglet s tat es e.g . 4sp2s 2 are insignificant. We divide the S^ z (q,u) DSF into its main contributing eigenstates in 
figures [5(d) | and [5 (e) These illustrations show clearly the difference in the momentum and energy localization of the 
2-spinon and 4-spinon signal. Moreover, although the weight distribution of the 4sp2oo is clearly broader in energy, 
the 2spls 2 signal is highly concentrated around q = ir and zero energy. The remarkable contribution of the 4-spinon 
clearly differentiates with the S zz (q, u) DSF at h = [50| 152] . 





(a) lpih 



(b) 2p2h 



(d) 2spls 2 
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(c) 2hls 2 
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0.14 

0.12 

0.1 

0.08 

0.06 

0.04 

0.02 



(e) 4sp2oo 



FIG. 5: Separated figures of the S^ z (q, uj} DSF with N = 400 cor resp onding to types of excitation in the labels. The 
magnetization is M — 150 in (a) (b) and (c) and M — 200 for (d) and (e) 



VIII. CONCLUSION 

Throughout this article we have described the con- 
struction of the SySj~ +1 and S Z S Z +1 form factors and we 
have used these results to formulate the adjacent spin op- 
erator correlation functions S h+ (g,Co>) and S^ z (q,uS). 
The sum rules of each DSF have been given and we have 
numerically evaluated the S h+ (g,cj) and S 4z (q,u) 
DSF in an isotropic spin chain with 400 sites in several 
magnetic fields. Among other results, we find a large 
contribution of the 4-spinon states in the S 4z (q,u) DSF. 
Moreover, we have shown that in some regions of the 
momentum-energy space, the signal can be attributed to 
4-spinon excitations and clearly distinguished from the 2- 
spinon signal. There is presently, no experimental setup 
which allows one to measure such an observable but we 
can trust that, with the progresses in realizations of spin 
chains with Bose gases in optical lattices [53-55 , mea- 



surements of adjacent spin operator correlation functions 
should be possible. 

Besides, the size scaling of Sj~ Sj~ +1 and S Z S Z +1 FFs 
at the spectrum lower boundary could be compared to 
results of the low energy effective field theory [56] for 
w«l and to the predictions of the non-linear Luttinger 
Liquid method [57] for higher energies. 

More generally, the results presented can be extended 
in different manners. The FFs of other operators Sf Sf+i 
with a, b = z, +, — can be constructed. FFs including 
more next neighboring spin operators should in principle 
also be possible keeping in mind that every S z or Sj~ op- 
erator adds a factor N or iV 2 , respectively, to the total 
computation load. Moreover, the present formula for the 
FFs have been expressed for a general anisotropy and can 
easily be implemented for the evaluation of S 4z (q,u) or 
S h+ (<7,cj) in the massless and massive regime, respec- 
tively A < 1 and A > 1. 
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